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1. Introduction 

There has been considerable interest in formulating measures of entanglement and 
applying them to extended quantum systems with many degrees of freedom, such 
as quantum spin chains, especially close to a quantum phase transition (see, e.g., 
several other reports in these proceedings). One of these measures is entanglement 
entropy. 1 Suppose the whole system is in a pure quantum state |\&), with density 
matrix p = and an observer A measures a subset A of a complete set 

of commuting observables, while another observer B may measure the remainder. 
A's reduced density matrix is pa — Tr# p. The entanglement entropy is just the 
von Neumann entropy Sa = — Tr^ pa log pa associated with this reduced density 
matrix. 

Although there are other measures of entanglement, the entropy is most readily 
suited to analytic investigation. Lack of space prevents us from discussing other 
measures of entanglement and for referring properly to all the works about en- 
tanglement entropy. For a more complete list of references we refer, e.g., to our 
papers. 2,3 

In these proceedings we try to give a pedagogical and non-technical introduction 
to the Quantum Field Theory (QFT) approach to entanglement entropy. 2 Particular 
attention is devoted to the one space dimensional case, with a linear dispersion 
relation, that, exactly at the quantum phase transition, can be effectively described 
by a two-dimensional Conformal Field Theory (CFT). Again, for lack of space, we 
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will not discuss the interesting topic of the time evolution of entanglement entropy 3 
and we will limit ourselves to discussion of the ground state properties. 

The layout of the paper is as follows. In the next section, we give an introduction 
to the concepts of renormalization group and conformal invariance. In section 3 wc 
give a path integral formula for the entanglement entropy. In Section 4, we apply 
this formula to the calculation of the entanglement entropy in 1+1 dimensional 
CFT. In section 5 we relax the conformal constraint and study the more general 
case of a system close, but not exactly at the phase transition. Several omitted 
(technical) details, as well as some results in higher dimensions, can be found in 
the original paper. 2 

2. Renormalization Group and Conformal Transformations 

For someone who is not familiar with the modern theory of phase transitions, the 
first natural question arising after reading the Introduction is what is the connec- 
tion between the QFT and a many-body system undergoing a continuous phase 
transition. The bridge between these apparently disconnected subjects is given by 
the Renormalization Group (RG) theory. Let us consider a quantum model de- 
fined by an hamiltonian H(g), where g is a tunable experimental parameter, e.g. 
in the well-known Ising model in a transverse field it is the magnetic field in the 
transverse direction. 5 Such hamiltonian at g = g c undergoes a continuous phase 
transition (g c is called Quantum Critical Point, QCP). Close to the QCP, the cor- 
relation length, that is the only relevant scale for the long-distance physics, behaves 
like £ ~ \g — g c \ ~ v , diverging at the QCP. [y > is an example of critical exponent.) 
Thus, at the QCP, the system is scale invariant. Now, the universality hypothesis 
states that some physical properties (called universal), close to the phase transition, 
do not depend on microscopic details, but only on global properties, such as sym- 
metries and dimensionality. This hypothesis has an elegant explanation in terms 
of RG theory: Under RG transformations (that, roughly speaking, are practical 
implementations of scale transformations), different hamiltonians sharing the same 
universal characteristics How to the same fixed point. This completely determines 
the long-distance behavior. For example a model defined on the lattice is defined by 
an hamiltonian which is invariant under translations multiple of the lattice spacing. 
The resulting fixed point hamiltonian, instead is generically invariant under arbi- 
trary translations, allowing for the use of a continuum field theory. For the same 
reason, the critical point hamiltonian is usually invariant under general rotations. 

These transformations of rotations, translations, and scaling form a group. Let 
us exploit the consequences of this symmetry group. Consider the two-point function 
of a scalar observable (0(ri)</>(r2)}. By translational invariance it can be only a 
function of ri — r2, by rotational invariance can depend only upon the modulus of 
such vector, and for a scale transformation r — > br it must behaves like 



(</»(r 1 )0(r 2 ))=6 2A *(0(&r 1 )0(6r 2 )), 



(1) 
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where the exponent is called scaling dimension of the field <j). These three con- 
ditions can be true if and only if 

(0(r 1 )0(r 2 )) = |r 1 -r 2 |- 2A *, (2) 

apart a normalization constant we set equal to 1. 

It turns out that a fixed point hamiltonian that is invariant under translations, 
rotations, and scaling transformations has usually the symmetry of the larger con- 
formal group 6 defined as the set of transformations that do not change the angles 
between two arbitrary curves crossing each other in some point. The consequences 
under this further invariance in a two-dimensional cuclidcan space 3, are extraordi- 
nary. In fact, in two-dimensions, we can use complex coordinates z = x + iy and 
z = x — iy. It is straightforward to prove that all the analytic functions f(z) are 
conformal transformations (this is the reason such mapping are so useful to solve 
Laplace's equations). The resulting symmetry group is infinite dimensional, and 
we can calculate, in principle, everything in an analytic way. For example, under 
a transformation of the type z — > w = w(z), we can generalize Eq. (1) to a space 
dependent scale factor b(z) = w'(z) (the prime here and in the following denotes 
the derivative with respect to z), obtaining 

(4>{zi)4>{z 2 )) = (\w'(z 1 )w'(z 2 )\) 2A ^ ( p(w(z 1 )) ( p(w(z 2 ))) . (3) 

This equation relates the two-point function of a scalar field on the plane given by 
Eq. (2) to the one in any other geometry (cylinder, torus, . . . ). 

A crucial rule in CFT is played by the so called stress tensor T^, that can be 
defined as follows. Under an arbitrary transformation x M — > x^ + e M , the euclidean 
action (the hamiltonian for classical systems) changes as 

S -> S + SS, with 5S = J d 2 xT^d„e v . (4) 

In general for a CFT T^ v can be chosen to be symmetric and traceless. b 

One of the most intriguing results 6 of two-dimensional CFT is that the univer- 
sality class (of minimal unitary models) is characterized by a single quantity called 
central charge, that assumes only the discrete values 

c=l — -, with m = 3, 4, . . . , oo . (5) 

m(m + 1) 

For example the Ising universality class corresponds to c = 1/2, the free boson to 
c = 1, and the three state Potts model to c = 4/5. Cases with c different from the 
value allowed by Eq. (5) are also of physical interest (e.g. the case c = is the critical 



a In a one-dimensional quantum system with a linear dispersion relation E oc k (this to ensure that 
space and time scale in the same manner), the two euclidean coordinates (x, y) correspond to the 
space and to the imaginary time r, obtained by Wick rotation of the real time t, i.e. y = r = it. 
b Noticc that in complex coordinates the traceless condition reads T Z z = 0. The standard conven- 
tion is T(z) = —2-kTzz and T(z) = —2-kTzz, called holomorphic and antiholomorphic stress tensor 
respectively. 6 
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percolation), but they are always pathological in some sense. From the knowledge 
of the central change c (that actually appears directly in several important physical 
quantities, as the entanglement entropy 7 ), we can (in principle) determine all the 
critical properties of the model, such as the critical exponents. 

3. Path integral formula for the Entanglement Entropy 

Consider a lattice quantum theory in one space and one time dimension. The lattice 
spacing is a, and the lattice sites are labelled by a discrete variable x. The domain 
of x can be finite, i.e. some interval of length L, semi-infinite, or infinite. Time is 
considered to be continuous. A complete set of local commuting observables will be 
denoted by {(j>(x)}, and their eigenvalues and corresponding eigenstates by {4*{x)} 
and <S> x \{4>(x)}) respectively. For a bosonic lattice field theory, these will be the 
fundamental bosonic fields of the theory; for a spin model some particular compo- 
nent of the local spin. The dynamics of the theory is described by a time-evolution 
operator H. The density matrix p in a thermal state at inverse temperature (3 is 

Piwrmw}) = zipy^wx'r^e-^mx 1 ) 1 }) , (6) 

where Z{[3) = Tre _/3ff is the partition function. 

This may be expressed in the standard way as a (cuclidcan) path integral: 

P = Z~ 1 f [dcf>(x, t)] J] 6(<P(x, 0) - <P(xJ) [] SMx, (3) - cfiix")") e- 5 * , (7) 

J X X 

where Se = Jq LecIt, with Le the euclidean lagrangian. 

The normalization factor of the partition function ensures that Trp = 1, and is 
found by setting {(j){x)"} = {(j){x)'} and integrating over these variables. This has 
the effect of sewing together the edges along r = and r = (3 to form a cylinder of 
circumference (3. 

Let A be a subsystem consisting of the points x in the disjoint intervals 
(ui,vi), . . . , (mat, wat). An expression for the the reduced density matrix pA may 
be found by sewing together only those points x which are not in A. This will leave 
open cuts, one for each interval (uj,Vj), along the the line r = 0. See figure 1 for a 
pictorial representation of this. 

We may then compute Tr//^, for any positive integer n, by making n copies 
of the above, labelled by an integer k with 1 < k < n, and sewing them together 
cyclically along the the cuts so that 4*{x)' k = (j){x)'l, +1 (and (f>(x)' n = (f>(x)") for all 
x e A. Let us denote the path integral on this n-sheeted structure by Z n (A). Then 

^P n A = ^- (8) 

Since Trp^ = . A™, where {\j} are all the eigenvalues of pa (which lie in [0, f)), 
it follows that the left hand side is absolutely convergent and therefore analytic for 
all Re n > f . The derivative with respect to n therefore also exists and is analytic in 
the region. Moreover, if the entropy pa = — J2\ ^ 1°S ^ i s finite, the limit as n — > f + 



I 
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of the first derivative converges to this value. We conclude that the right hand side 
of (8) has a unique analytic continuation to Ren > 1 and that its first derivative 
at n = 1 gives the required entropy: 

q ,. d d Z n (A) ... 
S A = - hm — Trp A = - hm — — . (9) 

So far, everything has been in the discrete space domain. We showed 2 that the 
continuum limit can be taken safely, since most of the UV divergences of the QFT 
cancel in the ratio (8). 



4. Entanglement entropy in 1+1-dimensional CFT. 

Now specialize the discussion of the previous section to the case when the field 
theory is relativistic and massless, i.e. a CFT, with central charge c. 



4.1. Single interval 

We first consider the case of a single interval of length £ in an infinitely long Id quan- 
tum system, at zero temperature. The ratio (8) is given by (0|0}7j. n , i- e - the vacuum 
expectation value in the n-sheeted surface. Within CFT to obtain this expectation 
value is enough to know how it transforms under a general conformal transforma- 
tion. This is formally given by (T(w))-jz n , where T(w) is the (holomorphic) stress 
tensor. The reason for that will be clear (we hope) in the following. 

To obtain (T(w))fi n , we need to map the n-sheeted surface onto a geometry 
where the mean value of the stress tensor is known, and then use the transformation 
law 6 

ill I _ 3 _//2 

T(w) = (z") 2 T(z) + ± Z - / . (10) 
12 z' 

The mapping we need is w — > z(w) = ((w — u)/(w — u)) 1 ^". This maps the 
whole of the n-sheeted surface lZ n to the z-plane C, where by translational and 
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rotational invariance (T(z))c = 0. In fact, the conformal transformation w — > ( = 
(w — u)/(w — v) maps the branch points to (0, oo). This is then uniformised by the 
mapping £ — > z = ( 1 ^ n . 

Thus, taking the expectation value of (10) and calculating the three derivatives, 
we find 

c(l-(l/n) 2 ) {v-uf 
(T(w))n n - ^ {w _ u) 2 {w _ v) 2- (H) 

Now the way of working should be clear. Compare this with the standard form of 
the correlator of T with two (primary) operators $ n (it) an d $_ n (u) which have the 
same complex scaling dimensions A„ = A„ = (c/24)(l — (1/n) 2 ) (conformal Ward 
identity): 6 

(T(w)$ n (u)<f>- n (v)} C = 7 ^7 A " x 2A _ 2 ,- rr^- , (12) 

(w — uY(w — v) z (v — uy^ n z (v~u) zlx ™ 

where $±„ are normalized so that ($ n (u)$_„(w))c = \v — u| _4A ". In writing the 
above, we are assuming that w is a complex coordinate on a single sheet C, which 
is now decoupled from the others. We have therefore shown that 

J[dd>]T(w)e- s ^) (mW$-„H>c m , 
[i(w ))nn - m]e s B (K n ) ($„(«)$_„(t;))c ' 1 j 

The insertion of T(w) on each sheet is given by (12). Since this is to be inserted 
on all the sheets, the right hand side gets multiplied by a factor n. Since the Ward 
identity (12) determines all the properties under conformal transformations, we 
conclude that Z n (A)/Z n cx Tr p n A behaves (apart from a possible overall constant) 
under scale and conformal transformations identically to the nth power of two- 
point function of a primary operator with A„ = A„ = (c/24)(l — (1/n) 2 ). In 
particular, this means that 

/,,_ „\ -(c/6)(n-l/n) 

Trpl = c„ , (14) 



where the exponent is just 4nA„. The power of a has been inserted so as the make 
the final result dimensionless, as it should be. The constants c„ are not determined 
by this method. However c\ must be unity. Differentiating with respect to n and 
setting n = 1, we recover the result of Holzhey et al. 7 

SA^log^+c;. (15) 

Notice that the constant c[ is not universal. 

The fact that Trp^ transforms under conformal transformations as a 2-point 
function of primary operators <fr± n means that it can be simply computed in other 
geometries, obtained by a conformal mapping z — > w — w(z), using equation (3). 
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4.2. Finite Temperature 

The transformation w — ► w' = {(3 /2tt) log w maps each sheet in the w-plane into 
an infinitely long cylinder of circumference (3. The sheets are now sewn together 
along a branch cut joining the images of the points u and v. By arranging this to 
lie parallel to the axis of the cylinder, we get an expression for Trp^| in a thermal 
mixed state at finite temperature /3 _1 . This leads to the result for the entropy 

^(/3)~^logfA sinh ^ +c ;. (16) 

For I <C (3 we find S A ~ ( c /3) log(£/a) as before, while, in the opposite limit i » /3, 
Sa ~ (ttc/3)(£/ (3). In this limit, the von Neumann entropy is extensive, and its 
density agrees with that of the Gibbs entropy of an isolated system of length I. 

Notice that the entropy given by Eq. (16) is a measure of quantum entanglement 
only for (3 — > oo. In the opposite limit (3 <C i it is just a measure of the classical 
entropy. Eq. (16) tells us how the crossover between these two different objects is 
realized increasing the temperature. 

4.3. Finite Systems 

Orienting the branch cut perpendicular to the axis of the cylinder (with the re- 
placement (3 — > L) corresponds to the entropy of a subsystem of length I in a finite 
system of length L, with periodic boundary conditions. This gives 

S A ~ ^logf — sin^-) +c[. (17) 
o \na L J 

Note that this expression is symmetric under i — > L— I. It is maximal when t = L/2. 

4.4. Finite system with a boundary. 

Next consider the case when the Id system is a semi- infinite line, say [0, oo), and the 
subsystem A is the finite interval [0, £). At x = a boundary condition is imposed, 
that must be conformal invariant as well. The n-sheeted surface then consists of n 
copies of the half-plane x > 0, sewn together along < x < £, r = 0. Once again, 
we work initially at zero temperature. It is convenient to use the complex variable 
w = t + ix. The uniformising transformation is now z = {{w — i£)/(w + it)^^ n , 
which maps the whole n-sheeted surface to the unit disc \z\ < 1. In this geometry, 
(T(z)) = by rotational invariance, so that, using (10), we find 

A (2f) 2 

( T H)n n = t . "\ ' W2 , (18) 

(w - ity(w + ity 

where A„ = (c/24)(l — n~ 2 ) as before. (18) has the same form as (T(w)$ n (i£)), 
which follows from the Ward identities of boundary CFT, 6 with the normalization 

<$ n (i*)) = (2^)- A ". 
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The analysis then proceeds in analogy with the previous case. We find 

Trpl^^e/a)^ 12 ^- 1 /^, (19) 

so that S A ~ (c/6) log(2£/a) + c[. 

Once again, this result can be conformally transformed into a number of other 
cases. At finite temperature (i~ x we find 

S A (0) ~ (c/6) log (09/tto) sinh(27rf//3)) + . (20) 

By taking the limit when <>/?we find the same extensive entropy as before. 

For a completely finite Id system, of length L (with two conformal boundary 
conditions at x = and x = L), at zero temperature, divided into two pieces of 
lengths I and L — £, we similarly find 

5 A = (c/6) log ((2L/na) sm^t/L)) + ci . (21) 
4.5. General case 

For the general case, when A consists of several disjoint intervals (uk,Vk), the 
uniformising transformation is rather complicated, 2 but the method is the same. 
After a quite long algebra we found 

c I V^i V k-Uj ^ u k -Uj ^ Vfc-^l . , /orA 
S *=3 \l^ lo S^^-2^ lo &^^-2^ l °S^^ +Nc,. (22) 

\j<fc 3<k j<k ] 

A similar expression holds in the case of a boundary, with half of the Wi corre- 
sponding to the image points (this is relevant in considering the time evolution, 3 
where the boundary, in the time axis, is the initial condition). 

5. Non-critical 1+1-dimensional models 

So far, we considered the case of a system that is exactly at the quantum critical 
point. The natural question arising is what happens when a system is close to the 
phase transition, in the region where the correlation length £ is large, but finite. 
In this case the system is still effectively described by a QFT, that is massive, the 
mass being the inverse of the correlation length. 

When the subset A is the whole negative real axis (and B the positive one), we 
showed that the entanglement entropy is Sa = (c/6) log£/a. The argument we gave 2 
is quite technical and we not repeat it here. However with this formula at hand, 
we can easily generalize it to more complicated situations, i.e., when A consists of 
several disjoint intervals, having A boundary points with its complement. In fact, 
when the correlation length is large, but smaller than the length of each interval, 
the disjoint pieces must not interact (this is nothing but cluster decomposition in 
ordinary QFT), resulting in 

S A = A\ log i . (23) 
6 a 



I 



I 
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Fig. 2. Entanglement entropy for the ID Ising chain as function of A. 



When the correlation length exceeds the lengths of all the intervals we crosses over 
to the CFT results. In the intermediate regime (i.e. £ of the order of some interval 
length), we expect to see very complicated crossover scaling forms, that are still 
universal, but probably depending on the model in a very complicated way, not 
only through the central charge. We emphasize once again that this result is true 
only when £ ^> a. 

Equation (23) is completely general, and, at the time of publication of our work, 
there was no evidence for its validity (apart from numerical results 8 suggesting the 
saturation of Sa for finite £). For this reason we checked this formula in several 
models in the geometry with A = 1. We considered the massive Gaussian QFT, 
reproducing Eq. (23) with c = 1. In this particular case we also calculated 2 the 
entanglement entropy in a finite system of total length 2L divided into two equal 
parts, showing explicitly the crossover between the CFT and the massive regimes. 

Then, relaxing even the condition £ 3> a, we solved exactly the Quantum Ising 
chain in a transverse field (by means of Corner Transfer Matrix method 9 ) with 
hamiltonian 

n n 

that displays a quantum phase transition for A = 1. The final result we obtained 
is given by complicated (but rapidly converging) infinite sums. 2 The resulting en- 
tanglement entropy as function of A is shown in Figure 2, in particular it diverges 
at A = 1, as it should. For £ = |A — ;» 1, it behaves as Sa ~ 1/12 log£, in 
agreement with Eq. (23) with c = 1/2. We also considered the XXZ model, 2 in the 
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vicinity of the isotropic point, confirming Eq. (23) for c = 1. 

Later, Peschel 10 was able to write the infinite sums in terms of elliptic integrals 
and also to generalize the result to the XY chain in a transverse field. Then a direct 
comparison with the explicit formula for A = 2 (obtained in the meantime by Its et 
al. 11 using a completely different approach) was possible, showing also the validity 
of the cluster decomposition. 
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